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U-Pb  zircon  geochronology  is  hampered  by  problems  acquiring  meaningful  geologic  ages  on  zoned  grains 
that  retain  isotope  signatures  from  multiple  growth  or  thermal  events.  We  present  a  new  method  using  laser 
ablation-multicollector-inductively  coupled  plasma-mass  spectrometry  to  overcome  complications  asso¬ 
ciated  with  intricately  zoned  zircon  crystals  through  in  situ  sampling  of  zircon  volumes  as  small  as  12-14  pm 
in  diameter  by  4-5  pm  in  depth  (<3  ng  of  zircon).  Using  Channeltron  multipliers  to  monitor  Pb  intensities  in 
conjunction  with  a  total  ion  counting  method  and  errors  calculated  as  function  of  the  number  of  counts,  the 
small-volume  technique  reproduced  published  ages  on  eight  Mesoproterozoic-Cretaceous  secondary  zircon 
standards  precise  and  accurate  within  2%,  and  an  age  ~  1  Ma  too  young  on  a  Oligocene-aged  grain.  Two  initial 
applications  of  the  small-volume  technique  —  the  detrital  zircon  provenance  of  fine-grained  mudstones  and 
shales  and  the  creation  of  zircon  U-Pb  age  maps  to  investigate  the  detrital  and  metamorphic  history  of  a 
granulite-facies  paragneiss  —  demonstrate  the  utility  of  this  technique  to  a  variety  of  geologic  problems  and 
confirm  the  viability  of  laser  ablation-multicollector-inductively  coupled  plasma-mass  spectrometry  as  a  tool 
for  high  spatial  resolution  U-Pb  geochronology. 

Published  by  Elsevier  B.V. 


1.  Introduction 

U-Pb  zircon  geochronology  has  become  one  of  the  most  powerful 
techniques  for  obtaining  absolute  ages  on  a  variety  of  high- 
temperature  igneous  and  metamorphic  events,  as  well  as  placing 
constraints  on  the  timing  and  provenance  of  sedimentary  processes. 
Zircon,  in  particular,  has  proven  especially  well  suited  for  U-Pb 
geochronology  because  of  its  relative  abundance  in  a  wide  variety  of 
rock  types,  its  tendency  to  include  U  and  exclude  Pb  from  its  crystal 
structure,  and  its  resistance  to  physical  and  chemical  weathering.  In 
addition,  zircon's  high  closure  temperature  to  the  diffusion  of  Pb 
allows  single  crystals  to  retain  chemically  distinct  zones  that  record 
multiple  growth  or  high-temperature  Pb-loss  events  (Cherniak  and 
Watson,  2003).  While  these  complexly  zoned  crystals  with  inherited 
cores  and  younger  rims  —  the  product  of  either  renewed  crystal 
growth  or  Pb-loss  —  contain  tremendous  potential  for  unraveling 
complicated  tectonic  histories,  these  grains  are  analytically  challen¬ 
ging  because  analysis  of  multiple  domains  can  yields  geologically 
meaningless  ages  that  represent  a  mixture  of  multiple  thermal  events. 
Although  the  decay  scheme  of  U,  which  includes  two  radiogenic 
isotopes  with  well-known  decay  constants,  facilitates  checks  for  the 
presence  of  mixed  ages  and  a  method  for  the  deconvolution  of 
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multiple  ages  (Wetherill,  1956),  the  confidence  in  deconvolved  older 
and  younger  ages  is  often  poor.  New  research  attempting  to  eliminate 
these  problems  with  complexly  zoned  zircons  takes  advantage  of 
recent  improvements  in  the  sensitivity  and  analytical  capabilities  of 
mass  spectrometers  by  analyzing  increasingly  smaller  samples, 
obtained  by  breaking  individual  crystals  (e.g.,  Schmitz  and  Bowring, 
2000)  or  using  in  situ  techniques  including  ion  beams  (e.g.,  Vavra 
et  al,  1996;  Ireland  and  Williams,  2003)  and  laser  ablation  (e.g.,  Eryer 
et  al.,  1993;  Rosier  and  Sylvester,  2003,  Gehrels  et  al.,  2008),  to  isolate 
and  selectively  sample  parts  of  zircons  that  are  characterized  by 
single-age  domains.  Here,  we  present  a  new  technique  using  laser 
ablation-multicollector-inductively  coupled  plasma-mass  spectrome¬ 
try  (LA-MC-ICP-MS)  at  the  Arizona  LaserChron  Center  (University  of 
Arizona)  to  sample  and  analyze  U-Pb  isotopes  from  zircon  volumes  as 
small  as  12-14  pm  in  diameter  by  4-5  pm  in  depth  (<3  ng),  and 
capable  of  rapidly  producing  ages  with  precision  and  accuracy  <2%  (at 
2a). 

2.  Small-volume  U-Pb  method 

2.1.  Sample  preparation 

Zircons  are  separated  from  bulk  rock  hand  samples  using  standard 
rock  pulverization  followed  by  density  and  magnetic  separation 
techniques.  Separated  zircons  are  then  mounted  with  an  external 
zircon  standard  in  epoxy  discs  one  inch  in  diameter;  mounting 
unknown  zircons  and  standards  close  to  each  other  and  in  the  inner  1  /2 
inch  of  the  mount  minimizes  U-Pb  fractionation  related  to  spatial 
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Fig.l.  Sketch  diagram  of  the  New  Wave  Instruments  laser-GV  Instruments  Isoprobe  at 
the  Arizona  LaserChron  Center  that  illustrates  the  collector  configuration  used  during 
small-voIume  U-Pb  geochronology. 


variation  of  carrier  gas  over  the  sample  surface  (Gehrels  et  al,  2008). 
Sample  mounts  are  then  sanded  and  polished  ~  1  /3  of  the  way  through 
unknown  grains,  carbon  coated,  and  imaged  on  a  SEM  fitted  with  a 
cathodoluminescence  (CL)  detector  to  characterize  chemical  zoning, 
qualitatively  assess  U  concentration,  and  identify  core-rim  relation¬ 
ships  within  unknown  grains  (Hanchar  and  Miller,  1993;  Corfu  et  al, 
2003;  Nasdala  et  al,  2003).  After  CL  imaging,  the  carbon  coat  and  any 
associated  common  Pb  is  removed  from  the  sample  surface  by  lightly 
repolishing  (removing  the  uppermost  ~l-2  pm)  and  ultrasonically 
washing  the  mount  in  a  dilute  acid  (1%  HN03  +  1%  HCl)  solution. 

2.2.  Data  acquisition:  laser  and  ICP-MS  setup 

U-Pb  analyses  are  performed  at  the  Arizona  LaserChron  Center  at  the 
University  of  Arizona  using  a  GV  Instruments  Isoprobe  equipped  with  the 
S-option  interface  and  coupled  with  a  New  Wave  Instruments  193  nm 
ArF  excimer  laser  (Fig.  1 ).  During  the  30-second  data  acquisition  routine, 
zircons  embedded  in  the  sample  mount  are  micro-sampled  in  situ  using 
a  laser  operating  with  the  beam  diameter  set  to  10  pm  and  a  constant 
fluence  of  ~4  J/cm^  (Gehrels  et  al.,  2008).  Laser  pits  imaged  and 
measured  using  a  MicroXAM  3D  profiler  at  the  University  of  California, 
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Fig.  3.  The  relationship  between  average  counts  per  second  and  total  integrated  counts 
^°®Pb  on  standard  Sri  Lanka  Zircon  using  24  laser  pulses  with  variable  A)  carrier  gas  flow 
and  B)  laser  hit  pulse  rate. 


Los  Angeles  indicate  approximately  cylindrical  laser  pits  with  openings 
of  12-14  pm  in  diameter  tapering  to  bases  8-10  pm  in  diameter  with 
<1  pm  of  floor  relief  (Fig.  2).  Adjusting  the  total  number  of  laser  pulses 
from  32  to  40  resulted  in  pit  depths  of  4.0  pm  and  5.0  pm,  respectively, 
yielding  a  constant  ablation  rate  of  0.125  pm/laser  pulse  (0.5  pm/second. 
Fig.  2)  and  a  sample  mass  of  <3  ng.  A  series  of  tests  were  run  on  a  single 
fragment  of  standard  Sri  Lanka  Zircon  (Gehrels  et  al.,  2008)  to  optimize 
the  laser  pulse  rate  and  the  rate  of  He-Ar  carrier  gas  (Gunther  and 
Heinrich,  1999)  flowing  over  the  sample  surface  (used  to  transport 
ablated  material  into  the  plasma).  Holding  the  laser  hit  rate  constant  at 
4  Hz,  increasing  carrier  gas  flow  from  0.2-0.36  L/min  yielded  a  linearly 


Fig.  2.  Images  and  cross-sectional  profiles  produced  using  the  UCLA  MicroXAM  3D  profiler  illustrating  the  shape  and  size  of  the  laser  pits  produced  in  zircon  using  A)  32  laser  pulses 
and  B)  40  laser  pulses. 


Please  cite  this  article  as:  Johnston,  S.,  et  al..  Small-volume  U-Pb  zircon  geochronology  by  laser  ablation-multicollector-ICP-MS,  Chemical 
Geology  (2008),  doi:10.1016/j.chemgeo.2008.11.004 


ARTICLE  IN  PRESS 


S.  Johnston  et  al  /  Chemical  Geology  xxx  (2008)  xxx-xxx 


3 


Fig.  4.  Graph  illustrating  the  relative  timing  of  the  background  measurement,  laser 
engagement,  the  total  count  integration  of  sample  signal,  and  the  corresponding  isotope 
intensities  over  the  course  of  a  typical  analysis.  Grey  boxes  indicate  the  duration  of  the 
laser  engagement  for  32  pulses  (dark  grey)  and  40  pulses  (light  grey).  Average 
backgrounds  and  peak  signal  intensities  are  given  in  counts  per  second 
(100000  cps  =  l  mV)  from  an  analysis  of  Sri  Lanka  zircon  during  the  February  standard 
session  (40  laser  pulses  at  4  Hz). 


increasing  relationship  between  average  counts  per  second  and 
integrated  total  counts  over  the  entire  30  s  analysis  (Fig.  3A,  see  data 
reduction  and  description  of  total  counts  technique  below).  This 
relationship  may  be  related  to  either  increased  pit-extraction 
efficiency  or  greater  elemental  evaporation  in  the  carrier  gas  at 
higher  carrier  gas  flows.  In  contrast,  while  holding  the  carrier  gas 
flow  constant  at  0.24  L/min,  increasing  laser  pulse  rate  from  2  to 
8  Hz  yielded  higher  average  counts  per  second,  but  total  counts 
remained  constant,  although  less  scatter  was  observed  at  lower  hit 
rates  (Fig.  3B).  Under  normal  operating  conditions,  the  laser  is  set  to 
pulse  at  4  Hz  with  a  carrier  gas  flow  of  0.24  L/min  in  order  to 
minimize  scatter  and  generate  signal  intensities  best  suited  to  the 
Isoprobe  detectors  from  zircons  with  a  range  of  U-Pb  concentra¬ 
tions,  although  both  of  these  settings  can  be  adjusted  to  accom¬ 
modate  samples  with  unusually  high  or  low  Pb  isotope 
concentrations. 

U-Pb  isotopes  are  counted  on  the  Isoprobe  using  a  30-second 
acquisition  routine  that  measures  backgrounds  followed  by  sample 
peaks  in  one-second  integration  periods.  Background  intensities  are 
measured  with  the  laser  off  for  the  initial  10  s  of  the  analysis,  after 
which,  the  laser  is  fired  for  8-10  s,  the  sample  is  introduced  to  the 
plasma,  and  sample  signal  is  monitored  for  the  remainder  of  the 
analysis  routine.  The  sample  signal  arrives  at  the  detectors  4-5  s  after 
the  laser  is  fired  and  washes  out  several  seconds  before  the  end  of  the 
30-second  analysis  (Fig.  4).  The  Isoprobe  is  operated  with  the 
intermediate  gas  at  1.0  L/min,  the  Ar  coolant  gas  at  14.0  L/min,  an  Ar 
flow  rate  of  0.24  mL/min  through  the  hexapole  collision  cell,  and  with 
an  accelerating  voltage  of  ~6  kV  (Gehrels  et  al.,  2008).  To  thermally 
stabilize  the  plasma  during  analyses  and  thereby  increase  the 
sensitivity  and  the  stability  of  the  signal,  analyses  are  run  in  a  wet 
plasma  created  by  continuously  aspirating  MilliQ  water  into  the 
plasma  using  a  microconcentric  nebulizer  with  an  uptake  of  50  pL/min 
coupled  with  an  Ar  flow  rate  of  0.340  L/min  (Gehrels  et  al.,  2008).  At 
these  operating  conditions,  approximate  background  levels  (mea¬ 
sured  by  scanning  over  the  peaks  using  a  Channeltron  detector  and 
depending  on  the  particular  analysis  session)  are  approximately  300  cps 
204pb+204Hg,  1700  cps  2°®Pb,  2000  cps  ^^^Pb,  140  cps  ^^^Th  and  30  cps 
With  the  laser  activated  at  4  Hz  and  ablating  at  ~  0.5  pm/s,  peak  ion 


intensities  on  standard  Sri  Lanka  zircon  (563.5  Ma,  -518  ppm  U, 
-118  ppm  Th,  206pb/204p|3^^^gQQQ^  Gehrels  et  al.,  2008)  are  typically 
<  50  cps  above  background  on  2°4pb,  igo  000  cps  2o®pb,  15  000  cps  ^^^pb, 
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Fig.  5.  Results  from  Channeltron  linearity  experiments  comparing  Pb  signal  variation 
(blank  —  corrected)  on  the  channeltron  detectors  with  respect  to  measured  on 
faraday  detector  H6.  White  lines  and  grey  areas  indicate  the  weighted  average  of  all 
analyses  ±1  standard  error.  Dashed  ellipses  in  C  and  D  represent  uncorrected  values 
whereas  solid  ellipses  represent  corrected  values. 
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460000  cps  and  3500000  cps  this  translates  to  a  U 
sensitivity  of  -6800  cps/ppm  U  (Fig.  4). 

Low  Pb  isotope  yields  during  small-volume  U-Pb  analyses  require 
that  ^^'^Pb,  2°®Pb,  and  ^^^Pb  are  measured  on  low-side  Channeltron 
electron  multipliers,  L7-L5,  while  ^^^Th  and  are  measured  on 
high-side  faraday  detectors  H5  and  H6  with  10^^  ohm  resistors  (Fig.  1 ). 
Channeltron  detectors  are  run  in  pulse-counting  mode  and  are  set  to 
trip  at  -1x10®  cps  while  faraday  detectors  can  easily  accommodate 
1x10®  cps.  Channeltron  linearity  was  checked  using  the  detector 
geometry  above  and  a  dead  time  correction  of  zero  ns  for  detector  L7, 
and  20  ns  for  detectors  L6  and  L5,  to  analyze  a  radiogenically  enriched 
U-Pb  solution  (^^^pb;  206p|3.  207p|3.  238y  ^^-^23:  1600)  system¬ 
atically  diluted  to  produce  the  range  of  Pb  signal  intensities  typically 
observed  in  zircon  during  laser  analyses  (Electronic  Supplement,  Table 
ESI  A,  B).  Each  solution  was  analyzed  twice  using  a  counting  routine 
with  30  one-second  integrations  to  calculate  the  mean  and  standard 
deviation  on  each  peak.  With  respect  to  detectors  L7  (^°^Pb)  and 
L6  (^°®Pb),  were  linear  at  the  95%  confidence  level  to  1600  cps  and 
750000  cps,  respectively  (Eig.  5A,  B).  Although  there  was  a  14% 
decrease  in  the  ^°^Pb/^®®U  at  intensities  <4000  cps  and  a  2-5% 
increase  in  ^°^Pb/^®®U  at  intensities  <10000  cps,  corrections  were  not 
made  to  these  detectors  because  the  measurements  were  within  error 
of  the  mean.  In  contrast,  detector  L5  (^°^Pb)  records  a  - 1%  decrease  in 
signal  from  5000-25000  cps  followed  by  a  gradual  8%  increase  in 
signal  from  25000-400000  cps  (Eig.  5C).  This  dip  in  signal  intensity 
on  the  L5  detector  is  corrected  offline  using  a  visual-fit  model  that, 
with  respect  to  yields  linear  relationships  at  the  95%  confidence 
level  to  400000  cps  (Eig.  5C).  As  the  Channeltron  detectors  age  with 
use  and  time,  it  will  be  essential  to  repeat  these  linearity  experiments 
periodically  to  check  linearity  and  recalibrate  the  correction  models,  if 
necessary. 


requires  that  sample  errors  are  calculated  as  a  function  of  the  number 
of  counts  rather  than  the  more  traditional  method  of  calculating  the 
error  on  a  series  of  isotope  ratios.  This  relationship  between  count  rate 
and  error  can  be  investigated  using  the  linearity  tests,  which  yield 
analytical  errors  derived  over  a  broad  range  of  concentrations  and  on 
all  detectors,  and  suggest  systematic  but  different  error-count  rate 
trends  for  Channeltron  and  faraday  detectors.  In  contrast  to  a  Poisson 
distribution  with  error  equal  to  root  n  —  commonly  assumed  for 
counts  of  large  numbers  —  the  linearity  tests  indicate  that  Channel¬ 
tron  errors  are  precise  to  -  2%  before  increasing  exponentially  at  count 
rates  <10000  cps,  whereas  faraday  errors  are  precise  to  -1.5%  before 
exponentially  increasing  at  count  rates  <350  000  cps  (Pig.  6,  Electronic 
Supplement,  Table  ESIA).  As  such,  the  analytical  error  on  each  one- 
second  integration,  Os,  is  empirically  calibrated  from  the  linearity 
experiments  and  proportional  to  on  channeltron  and  faraday 
detectors,  respectively  (normalized  by  a  factor  of  1.02,  the  student's  t 
critical  value  for  two-tailed  confidence  at  68.27%  (la)  and  29  degrees 
of  freedom): 


CTs.Channeltron  =  1  -02  ((n^  +  200)°'"2  +  0.02(ns  +  200))  (2) 

and, 

CTs, faraday  =  102((ns  +  100000)°=  +  0.015(ns  +  100000))  (3) 

The  analytical  error  on  the  total  counts,  is  then  calculated  by 
adding  Os  to  a^g  in  quadrature  for  integrations  12-30: 


Channeltron 


I 

t=  12 


O': 


s, Channeltron 


+  a] 


bg 


(4) 


2.3.  Data  reduction:  total  count  technique,  error  propagation  and 
fractionation 

U-Pb  isotope  data  collected  on  the  Isoprobe  are  reduced  offline 
using  an  Excel  macro  and  spreadsheet  that  calculates  backgrounds 
and  performs  a  total  count  analysis  on  the  cumulative  sample  signal. 
Signals  collected  on  the  backgrounds  during  integrations  3-10  are 
used  to  calculate  average  background  intensities,  n^g,  and  correspond¬ 
ing  standard  deviation,  a^g  (normalized  to  the  Student's  t-distribution 
for  7  degrees  of  freedom),  on  each  of  the  measured  isotopes.  Por  ^®^Th 
and  background  levels  are  much  less  than  the  detection  limits  of 
the  faraday  detectors  on  which  they  are  measured,  and  n^g  and  a^g 
represent  a  measurement  of  the  electronic  baseline  for  each  individual 
detector.  In  contrast  to  the  relatively  constant  background  signals, 
signal  intensities  produced  on  the  sample  cannot  be  used  to  calculate 
meaningful  average  signal  intensities  due  to  a  number  of  factors 
including:  1)  noise  inherent  to  sample  introduction  through  laser 
ablation,  2)  decreasing  down-pit  extraction  efficiency,  3)  low  signals 
that  require  longer  counting  integration  periods  of  1  s,  and  4)  the  short 
8-10  s  duration  of  the  sample  signal.  To  avoid  this  noise,  signal 
intensities  on  the  sample  are  reported  as  a  total  count  integration,  N, 
of  the  counts  on  the  sample,  Hs,  minus  n^g,  for  analysis  integrations 
12-30  (Pig.  4): 


N  = 


t  =  30 


I  (ns-ntg) 

t  =  12 


(1) 


In  addition  to  eliminating  errors  associated  with  laser  signal  noise, 
another  advantage  of  the  total  counts  technique  is  that  it  eliminates 
the  need  to  correct  for  U-Pb  fractionation  associated  with  decreasing 
and  variable  down-pit  extraction  efficiency. 

While  the  total  count  method  successfully  avoids  errors  associated 
with  laser  signal  instability  and  down-pit  extraction  efficiency,  it 


and. 


faraday  “ 


(5) 


Although  this  method  for  error  propagation  cannot  detect  in-run 
variability  in  instrument  noise,  it  does  produce  a  robust  estimate  of 


Fig.  6.  Relationship  between  average  counts  per  second  (not  blank  —  corrected)  and  la 
error  from  sets  of  30  one-second  integrations  measured  on  Channeltron  and  faraday 
detectors  throughout  the  linearity  experiments.  Also  plotted  are  visual-fit  error 
calibration  curves  used  to  calculate  errors  during  laser  acquisition;  an  additional 
error  curve  based  on  root  n  is  plotted  for  reference. 
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Fig.  7.  Fractionation  corrected  A)  ^°®Pb/^‘’’Pb  and  B)  ^°®Pb/^^*U  ratios  measured  on 
Sri  Lanka  Zircon  standard  throughout  the  course  of  the  February  analysis  session 
using  40  laser  pulses  shown  with  measurement  errors  only.  Additional  required 
external  error  of  1.1%  per  analysis  on  206p|^y238y  interpreted  as  the  result  of  matrix 
affects  associated  with  ablation  characteristics,  depth-dependent  mass  fractiona¬ 
tion  combined  with  minor  variations  in  pit  shape,  or  subtle  variations  in  trace 
element  concentration. 


the  analytical  uncertainty  based  on  the  integrated  concentration  of  U 
and  Pb  isotopes  in  the  sample. 

Isotope  fractionation  —  the  cumulative  result  of  mass  fractiona¬ 
tion  due  to  variable  pit-extraction  and  ionization  efficiency, 
collector  gain  and  instrument  drift  —  is  corrected  using  fragments 
of  Sri  Lanka  Zircon  (563.5 ±3.2  Ma,  Gehrels  et  al.,  2008)  as  a 
standard  that  is  placed  in  the  sample  mount  adjacent  to  unknown 
grains  and  analyzed  every  4-5  analyses  as  an  external  standard. 
Measured  206p|^^207p|^  206p|^y238y  j-^^jos  on  the  standard  are 

compared  to  the  known  values  and  used  to  generate  time-sensitive 
fractionation  curves  from  which  running  averages  are  calculated 
and  applied  to  ratios  from  unknown  grains  (Gehrels  et  al.,  2008). 
Because  the  total  count  integration  method  used  in  the  small- 
volume  technique  does  not  make  a  correction  for  down-pit 
variations,  U-Pb  ratios  determined  by  this  method  are  more 
sensitive  to  down-pit  variations  in  extraction  efficiency  than  ratios 
determined  by  the  zero-time  intercept  method  (e.g.,  Gehrels  et  al, 
2008),  and  the  assumption  that  standards  and  unknowns  share 


similar  ablation  characteristics  is  particularly  important.  The 
propagation  of  errors  through  isotope  ratios  was  investigated 
using  fractionation-corrected  standard  analyses  from  a  typical 
analysis  session  (26  analyses  on  the  standard  over  the  course  of 
the  entire  session).  Adding  percent  errors  on  0^,206  and  0^,207 
(calculated  as  shown  above)  in  quadrature  yields  a  single  popula¬ 
tion  with  a  MSWD  <1  for  the  standards 

analyzed  throughout  the  session  (Fig.  7A).  In  contrast,  the  ^°®Pb/ 
ratios  on  standards  from  the  same  analysis  session  with  errors 
propagated  from  0^,206  and  0^,238  alone,  yields  a  large  MSWD  that 
requires  an  additional  1.1%  error  (la)  for  each  point  (Fig.  7B, 
external  error  after  Ludwig,  2003).  This  data  implies  that  a  source  of 
error  in  addition  to  the  measured  errors  on  ^ospj^  238y  _ 
related  to  either  matrix  affects  associated  with  ablation  character¬ 
istics,  depth-dependent  mass  fractionation  combined  with  minor 
variations  in  pit  shape,  or  subtle  variations  in  trace  element 
concentration  (Black  et  al.,  2004)  —  is  needed  to  describe  the 
206pb/238u  ratio.  As  such,  and  following  similar  patterns  on  Sri 
Lanka  Zircon  standard  ratios  from  other  analysis  sessions,  analytical 
errors  on  the  206p|^y207p|^  analysis  are  propagated 

from  af^^206  and  (Jf^^207  alone,  whereas  analytical  errors  on  the  ^°®Pb/ 
ratio  are  calculated  by  adding  the  measured  percent  errors  on 
^°®Pb  and  in  quadrature  to  a  user-added  external  error  of  1%. 

^°^Pb  is  monitored  during  analyses,  although  low  common  Pb 
concentrations  in  zircon  in  combination  with  small-volume  analysis 
typically  yields  count  rates  <50  cps  above  background,  correspond¬ 
ingly  high  errors,  and  values  that  are  not  significantly  different  than 
zero.  During  the  course  of  this  study,  ^°^Pb  was  used  primarily  to 
check  against  the  accidental  sampling  of  Pb-rich  inclusions  or  for 
surface  contamination,  and  a  common  Pb  correction  was  not  applied 
to  any  analyses. 

Approximate  U  and  Th  concentrations  are  calibrated  by  compar¬ 
ing  concentrations  in  the  Sri  Lanka  Zircon  (measured  by  isotope 
dilution  thermal  ionization  mass  spectrometry,  ID-TIMS: 
U  =  518  ppm,  Th  =  68  ppm,  Gehrels  et  al.,  2008)  to  the  average 
measured  signal  over  the  course  of  the  entire  session  analysis  to 
create  U  and  Th  concentration  correction  factors.  These  correction 
factors  are  then  used  to  generate  U  and  Th  concentrations  and 
determine  a  U/Th  ratio. 

Fractionation  corrected  ratios,  analytical  errors,  and  error  correla¬ 
tions  are  used  to  calculate  ages  and  assess  concordancy.  Typically,  for 
single  analysis  ages  (e.g.  ages  from  single  grains  in  a  detrital  sample), 
206pb/238u  3g05  analytical  errors  are  assigned  to  analyses  <1.2  Ga 
and  206p|^y207p|^  analytical  errors  are  assigned  to  analyses 

older  <1.2  Ga.  For  cogenetic  analyses  (e.g.,  zircon  analyses  interpreted 
to  have  shared  a  common  U-Pb  isotope  history),  reported  ages  are 
weighted  averages  of  either  206p|^y238y  ^g^^  samples  <1.2  Ga  or 
206pb/207pb  3gg5  fQj.  samples  older  <1.2  Ga,  whereas  all  errors  are 
quoted  as  the  2a  standard  error  on  the  mean  unless  stated  otherwise. 
More  complicated  samples  that  display  significant  inheritance  or  Pb- 
loss  are  typically  assessed  using  a  variety  of  factors  including  chemical 
zoning  with  respect  to  CL  imaging,  U  concentration  and  U/Th  ratio, 
and  ages  and  errors  are  determined  using  a  combination  concordia 
regressions,  weighted  averages  and  the  Tufflirc  algorithm  of  Ludwig 
(2003). 

Once  an  age  on  a  single  analysis  or  group  of  analyses  has  been 
determined,  systematic  errors  associated  with  the  age  must  be 
calculated.  2a  systematic  errors  present  in  the  small-volume 
technique  are  the  result  of  uncertainty  on  U  decay  constants  (0.16% 
for  238U  and  0.21%  for  235U,  Jaffey  et  al.,  1971;  Mattinson,  1987), 
uncertainty  on  the  age  of  the  external  standard  (0.57%,  Gehrels  et  al., 
2008),  and  the  average  uncertainty  on  the  running  average  deter¬ 
mined  through  fractionation  calibration  (usually  1.0-1.5%  depending 
on  the  analysis  session).  These  systematic  errors  are  added  in 
quadrature  to  analytical  errors,  and  yield  the  final  age  and  total 
error  for  the  sample. 
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3.  Secondary  zircon  standard  results:  precision  and  accuracy 


To  test  the  accuracy  and  precision  of  the  small-volume  technique,  a 
sample  mount  was  made  with  nine  secondary  zircon  standards  with 
ID-TIMS  ages  ranging  from  the  Mesoproterozoic-Miocene.  From 


oldest  to  youngest,  analyzed  secondary  standards  include:  FCl 
(1099.9  Ma,  Paces  and  Miller,  1993),  91,500  (1065  Ma,  Wiedenbeck 
et  al,  1995),  Peixe  (564  Ma,  Gehrels,  unpublished  data),  R33  (419.3  Ma, 
Black  et  al.,  2004),  Temora  (416.8  Ma,  Black  et  al.,  2004),  Plesocize 
(337.1  Ma,  Slama  et  al.,  2008),  49127  (136.6  Ma,  Mattinson  et  al.,  1986, 
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Fig.  8.  Concordia  plots  for  secondary  zircon  standards  analyzed  during  the  October  (open  ellipses)  and  February  (grey  ellipses)  small- volume  U-Pb  geochronology  sessions.  Dashed  ellipses 
denote  analyses  omitted  from  weighted  average  calculations.  Weighted  average  ages  and  total  2a  error  (^^^abbreviated  for  206p]^y238y  ^/^abbreviated  for  age,  ^/^abbreviated 

for  206p|^y207p|^  ^g^^  boxes  inset  within  the  concordia  diagrams;  ages  annotated  with  an  asterisk  indicate  ages  derived  from  data  with  an  MSWD  <  1.5  and  a  probability  of  fit 

<0.15. 
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Fig.  8  (continued). 


Kimbrough,  personal  communication),  Ecstall  (91.5  Ma,  Butler  et  al., 
2002),  and  Fish  Canyon  Tuff  (28.4  Ma  without  Th  correction,  Schmitz 
and  Bowring,  2001).  Secondary  standards  were  analyzed  in  two 
sessions  over  a  period  of  4  months:  using  7  analyses  per  sample  with 
32  laser  pulses  in  October,  2007  and  10  analyses  per  sample  with  40 
laser  pulses  in  February,  2008  (Fig.  8,  Electronic  Supplement,  Table 
ES2).  Average  ^°^Pb,  ^°®Pb,  and  ^°^Pb  backgrounds  from  individual 
analyses  measured  throughout  the  two  sessions,  respectively,  were 
410  cps,  2370  cps,  and  2580  during  the  October  session  and  200  cps, 
990  cps,  and  1500  cps  during  the  February  session;  electronic 
baselines  on  the  two  faraday  detectors  did  not  change  between  the 
two  analysis  sessions.  Background-corrected  signals  on  the  zircon 
samples  ranged  from  -1700-400000  cps  on  ~100- 

40  000  cps  on  ^°^Pb.  The  fractionation  correction  factor  ranged  from 
1.62-1.72  for  206p5/238y  1.18-1.26  for  slightly 

less  variation  in  the  Sri  Lanka  Standard  and  the  sliding  window  around 
the  fractionation  factor  during  the  October  session  yields  a  cumulative 
systematic  error  of  1.4%  as  opposed  to  1.5%  (2a  standard  error  on  the 
mean)  during  the  February  session. 

3.1.  Precision  and  reproducibility 

In  general,  analyses  from  both  sessions  yield  concordant 
analyses  that  overlap  within  error  of  each  other  (Fig.  8).  Out  of 


146  total  individual  analyses,  206p|^^238y  yielded  an  average 
2a  analytical  error  of  2.9%  and  an  average  2a  total  error  of  3.2%. 
207pb/235u  yielded  an  average  2a  error  of  7.7%  and  an  average 
2a  total  error  of  7.9%  (excluding  6  discarded  analyses  and  3 
additional  analyses  with  low  ^°^Pb).  Combining  multiple  analyses 
on  the  same  zircon  sample  and  from  the  same  analysis  session 
yields  weighted  mean  ages  with  precision  significantly  better  than 
individual  analyses  (Table  1).  Over  the  course  of  both  sessions,  6 
analyses  out  of  146  were  identified  as  outliers  using  either  clear 
discordance  or  2a  variation  from  the  mean  as  rejection  criterion, 
and  omitted  from  weighted  mean  calculations.  8  of  9  samples 
during  the  February  analysis  session  with  suitable  206p|^y238y 
MSWDs  <1.5  (high  MSWD  of  2.6)  and  6  of  8  samples  that  achieved 
the  same  level  of  reproducibility  (high  MSWD  of  2.0)  during  the 
October  session,  indicate  that  the  empirically  calibrated  user- 
added  external  error  of  1%  to  all  206p|^^238y  j-^^os  is  applicable  to  a 
wide  variety  of  zircons.  206p|^^207p|^  ratios  were  consistently 
reproducible  during  both  analysis  sessions  with  MSWDs  ranging 
from  0.6-1.4  in  February  and  from  0.5-1.4  in  October  (Table  1). 
Precision  on  weighted  mean  ages  was  slightly  better  during  the 
February  session  with  an  average  2a  total  error  of  1.8%  as  opposed 
to  1.9%  during  the  October  session.  It  is  possible  that  this  improved 
precision  and  reproducibility  observed  in  the  February  session  was 
the  result  of  better  sensitivity  and  increased  signal  to  background 
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Table  1 

U-Th  data  and  U-Pb  ages  on  secondary  zircon  standards  measured  by  LA-MC-ICP-MS  during  both  Ocotber  and  February  analysis  sessions 


Sample 

U 

(ppm) 

±la^ 

(ppm) 

U/Th 

±la^ 

206p|3/ 

238yb 

±2a 

(Ma)" 

±tot. 

(Ma)° 

MSWD 

207p|3/ 

235ub 

±2a 

(Ma)" 

±tot. 

(Ma)^ 

MSWD 

206p^l 

207pt,b 

±2a 

(Ma)" 

±tot. 

(Ma)° 

MSWD 

October  analysis  session:  32  laser  pulsesl4  pm  pits 
Peixe  207  92  8.6  2.58 

538.3 

5.9 

9.6 

1.40 

549.0 

10.0 

12.6 

0.86 

589 

45 

46 

0.55 

91,500 

109 

23 

2.7 

0.09 

1044.0 

11.0 

18.3 

0.92 

1043.0 

14.0 

20.2 

0.72 

1038 

35 

38 

0.53 

R33 

164 

49 

1.4 

0.18 

421.8 

5.3 

7.9 

0.77 

416.0 

12.0 

13.3 

1.06 

386 

73 

73 

0.94 

Temora 

131 

29 

2.8 

0.70 

419.9 

5.1 

7.8 

1.40 

422.0 

13.0 

14.3 

1.40 

445 

81 

81 

1.40 

49,127 

211 

60 

1.1 

0.09 

134.1 

1.8 

2.6 

1.50 

132.2 

9.9 

10.1 

1.13 

149 

180 

180 

1.11 

FCl 

529 

226 

2.1 

0.25 

1096.0 

11.0 

18.9 

2.00 

1102.4 

9.8 

18.3 

0.50 

1115 

20 

26 

0.56 

Fish  CT 

491 

364 

2.1 

0.73 

27.1 

0.5 

0.6 

2.00 

28.3 

3.5 

3.5 

0.40 

197 

290 

290 

0.47 

Plesovice 

685 

110 

11.6 

0.61 

333.3 

3.6 

5.9 

0.95 

329.3 

5.1 

6.9 

0.38 

303 

33 

33 

0.65 

February  analysis  session:  40  laser  pulses/5  pm  pits 
Peixe  173  77  9.0  2.37 

567.4 

5.1 

9.9 

1.15 

561.3 

6.8 

10.8 

0.83 

542 

27 

28 

0.61 

91,500 

85 

10 

2.9 

0.10 

1068.8 

9.2 

18.5 

1.30 

1056.2 

9.1 

18.3 

0.92 

1031 

21 

23 

0.83 

R33 

148 

62 

1.7 

0.45 

419.8 

3.9 

7.4 

0.75 

418.6 

6.4 

9.0 

0.44 

415 

35 

35 

0.98 

Temora 

148 

82 

3.2 

0.75 

413.2 

3.7 

7.2 

0.90 

408.0 

6.3 

8.8 

0.23 

386 

36 

36 

0.89 

49,127 

262 

70 

1.3 

0.34 

135.7 

1.2 

2.4 

0.94 

132.9 

3.5 

4.0 

1.05 

93 

62 

62 

1.40 

FCl 

419 

176 

2.2 

0.58 

1102.2 

9.1 

18.9 

1.40 

1108.0 

7.8 

18.4 

1.13 

1120 

15 

19 

0.66 

Fish  CT 

522 

277 

2.2 

0.95 

27.4 

0.3 

0.5 

1.20 

30.5 

1.5 

1.6 

0.70 

289 

110 

110 

0.98 

Ecstaii 

335 

144 

3.0 

0.77 

89.9 

1.5 

2.0 

2.60 

89.2 

2.7 

3.0 

0.44 

68 

71 

71 

0.74 

Plesovice 

832 

264 

10.3 

1.34 

335.7 

2.8 

5.8 

0.86 

335.8 

3.4 

6.1 

1.10 

338 

19 

19 

0.81 

^  U  concentration  and  U/Th  errors  are  1  standard  deviation. 

^  U-Pb  ages  are  weighted  averages  of  7  analyses  from  the  October  analysis  session  and  10  analyses  from  the  February  analysis  session.  U  decay  constants  and  composition  as 
follows:  238u=9.8485x10"^°  235u  =  1.55125x10“^°, 

^  2a  standard  error  on  the  weighted  average  age  accounting  for  analytical  errors  only. 

^  2a  standard  error  on  the  weighted  average  age  accounting  for  analytical  errors  and  systematic  errors;  accumulated  206p|^y238y  206pj^y207p|^  systematic  error,  respectively,  was 

1.4%  and  1.5%  during  the  October  analysis  session  and  1.5%  and  1.0%  during  the  February  analysis  session. 


levels.  However,  given  the  greater  reproducibility  of  the  Sri  Lanka 
Zircon  standard  and  the  corresponding  stability  of  the  fractiona¬ 
tion  factor  during  the  October  session,  it  is  likely  that  both  the 
greater  reproducibility  and  precision  observed  in  the  February 
session  are  simply  the  result  of  more  analyses  per  sample  during 
that  session. 

3.2.  Accuracy 

Ages  on  the  nine  secondary  zircon  standards  were  also  generally 
accurate  within  their  assigned  analytical  errors  with  respect  to 
published  ID-TIMS  ages.  Out  of  the  146  individual  analyses,  81%  of 
206pb/238u  3g0s  accurate  within  the  calculated  2a  total  error 
and  95%  of  analyses  accurate  within  5%,  whereas  91%  of  analyses 
yielded  ages  were  accurate  within  the  calculated  2atotal 

error  and  76%  of  analyses  were  accurate  within  5%.  Plotting  the 
207pb/235u  3g0  offset  from  the  ID-TlMS  age  for  all  analyses  against 
average  ^°^Pb  cps  during  laser  acquisition  (Fig.  9B)  does  not  reveal 
any  age  bias  with  respect  to  ^°^Pb  intensity  and  suggests  that  the 
linearity  correction  on  detector  L5  (^°^Pb)  is  robust  over  the  range  of 
Pb  concentrations  typically  observed  in  zircon.  In  contrast,  a  similar 
plot  illustrating  206p|^^238y  against  ^°®Pb  intensity 

displays  an  apparent  downward  trend  in  age  offsets  from  0  to  -4% 
over  206p|3  intensities  of  10  000-2000  cps,  and  indicates  that  the 
assumption  of  Channeltron  linearity  on  detector  L6  (^°®Pb)  may  be 
compromised  in  zircons  with  low  ^o^pj^  yields  (Fig.  9A).  Weighted 
mean  ages  were  accurate  within  2%  of  the  published  age  on  6  of  8 
samples  during  the  October  session  and  8  of  9  samples  during  the 
February  session;  7  of  9  samples  were  accurate  to  1%  during  the 
February  session  (Fig.  lOA,  Table  1).  With  respect  to  the  calculated 
total  2a  errors,  5  of  8  samples  during  the  October  session  and  8  of  9 
samples  during  the  February  session  were  accurate.  The  most 
problematic  sample  was  the  Fish  Canyon  Tuff  which  yielded  ages 
4.8%  and  3.6%  ( 1.3-1. 0  Ma)  too  young  in  the  October  and  February 
sessions,  respectively,  and  were  also  both  inaccurate  with  respect 
the  accepted  age  and  the  assigned  total  2a  errors.  These  inaccurate 
ages  are  likely  the  result  of  problems  with  Channeltron  linearity  at 


low  intensities.  In  particular,  the  increase  in  206p|^^238y  observed  in 
the  solution  linearity  experiments  suggests  that  2°®Pb  at  low 
intensities,  and  thus  the  ^oepj^  backgrounds,  were  likely  over¬ 
counted.  Although  overcounting  backgrounds  by  3-5%  (~50- 
85  cps)  has  little  effect  on  zircons  with  high  radiogenic  yields, 
overcounting  backgrounds  in  zircons  that  yield  low  2°®Pb  inten¬ 
sities  could  produce  erroneously  young  ages.  As  such,  the  young 
ages  observed  on  the  Oligocene-aged  Fish  Canyon  Tuff  (yielding 
-4700-6800  cps  ^°®Pb)  suggests  a  potential  limit  for  the  small- 
volume  U-Pb  technique  to  produce  precise  and  accurate  ages  (at 
better  than  2%)  in  zircons  with  ^ospj^  yields  <~io  000  cps  until 
Channeltron  linearity  on  detector  L6  is  better  calibrated  at  low 
intensities.  During  the  October  session,  young  and  inaccurate  ages 
on  Peixe  and  91500,  the  first  two  samples  analyzed,  are  likely  the 
result  of  rapid  instrument  drift  after  changing  the  sample  mount, 
and  associated  with  an  observed  exponentially  decaying  decrease 
in  backgrounds  of  19%  on  both  ^ospj^  207p|^  20 

analyses. 

3.3.  U-Th  concentrations 

U  and  Th  concentration  data  from  both  ID-TIMS  and  the 
small-volume  technique  presented  here  are  extremely  variable 
with  standard  deviations  up  to  50%  for  U  concentrations  and  up 
to  40%  for  U/Th  ratios.  Despite  this  large  intra-sample  variability, 
average  U  concentrations  and  U/Th  ratios  measured  via  ID-TIMS 
and  during  the  course  of  this  study  are  within  30%  of  each  other 
(Fig.  lOB,  C). 

4.  Applications 

4.1.  Detrital  zircon  provenance  of  mudstones  and  shales 

The  ability  to  rapidly  collect  precise  and  accurate  (<2%)  U-Pb 
ages  on  individual  zircon  grains  have  brought  LA-ICP-MS  to  the 
forefront  of  provenance  and  sediment  dispersal  studies  which 
investigate  detrital  zircon  age  populations  in  clastic  sedimentary 
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(average  cps) 


207 

Pb  (average  cps) 


Fig.  9.  A)  206p|3y238u  g)  207pgy235y  ID-TIMS  foF  ail  individual 

analyses  during  October  and  February  secondary  standard  analysis  sessions  plotted  as  a 
function  of  Pb  concentration.  Grey  boxes  indicate  error  windows  of  ±2  and  5%. 


rocks  (Gehrels  et  al.,  2006;  Gehrels  et  al,  2008).  In  these  studies, 
zircons  separated  from  a  sample  are  selected  randomly  and 
individually  dated  in  order  to  determine  a  detrital  age  probability 
distribution.  Age  distributions  can  then  be  used  to  determine  likely 
sediment  sources,  compared  to  distributions  from  sediments  in 
adjacent  basins  to  investigate  petrogenetic  history  and  sediment 
dispersal  patterns,  or  used  to  determine  maximum  depositional  age 
for  the  sampled  sediment.  However,  because  typical  LA-ICP-MS  labs 
use  laser  spots  25-50  pm  in  diameter  and  ablate  pits  12-20  pm 
deep  or  raster  across  the  grain  surface,  detrital  zircon  studies  using 
LA-ICP-MS  have  been  limited  to  coarse  siltstones  and  sandstones 
with  grain  sizes  <~30  pm  in  diameter.  The  small-volume  technique 
described  here  expands  the  range  of  grain  sizes  measurable  by  LA- 
ICP-MS  down  to  fine  silts  with  grain  sizes  <~  15  pm,  and  allows  for 
the  routine  characterization  detrital  zircon  age  populations  in 
siltstones  and  shales.  In  a  first  application  of  this  technique  to 
detrital  zircon  geochronology,  zircons  were  separated  from  the 
Cambrian  Bright  Angel  Shale  of  the  Grand  Canyon's  Tonto  Group. 
Detrital  zircons  separated  from  the  Bright  Angel  Shale  range  in  size 
from  20-50  pm  in  diameter  and  were  analyzed  using  40  laser 
pulses  with  a  spot  size  of  10  pm.  Out  of  105  analyses,  17  analyses 


were  discarded  due  to  either  extreme  discordance  or  high  errors, 
whereas  88  grains  yielded  approximately  concordant  ages  ranging 
from  1.0-2.4  Ga  (Fig.  11  A,  Electronic  Supplement,  Table  ES3A,  ES3B). 
The  age  probability  distribution  indicates  significant  peaks  at  1.0, 
1.4  and  1.7  Ga  with  only  two  Archean  grains;  including  discarded 
analyses  does  not  significantly  alter  the  distribution  (Eig.  11 B).  This 
age  distribution  contains  identical  age  peaks  to  the  Cambrian 
Tapeats  Sandstone  (Eig.  11 B,  Stewart  et  al.,  2001,  Gehrels  unpub¬ 
lished  data),  which  stratigraphically  underlies  the  Bright  Angel 
Shale,  and  supports  the  utility  of  this  technique  to  accurately 
characterize  detrital  zircon  age  populations. 

4.2.  U-Pb  age  mapping  of  complexly  zoned  zircon 

One  of  the  primary  goals  that  motivated  the  development  of  the 
small-volume  U-Pb  geochronologic  technique  was  to  investigate  U-Pb 
age  histories  of  zircon  populations  with  intricate  isotopic  zoning  at 
spatial  scales  <20  pm.  A  granulite-facies  pelitic  gneiss  collected  in 


A.  Age  offset  from  ID-TTMS  (%) 


Fig.  10.  Comparison  of  smali-voiume  U-Pb  geochronology  A)  age  dates,  B)  U  concentra¬ 
tions  and  C)  U/Th  values  with  published  ID-TlMS  data.  Grey  areas  in  B  and  C  represent  error 
windows  of  ±20%. 
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U-Pb  age  (Ma) 

Fig.  11.  A)  Concordia  plot  of  ail  analyzed  detritai  zircons  from  the  Bright  Angel  Shale. 
B)  Detritai  zircon  age  probability  distributions  from  the  Bright  Angel  Shale  compared 
to  the  distribution  from  the  immediately  underlying  Tapeats  Sandstone  U-Pb  age  is 
206pb/238u  for  ages  <1.2  Ga  and  206p|^y207pj^  foj.  ^ges  <1.2  Ga.  Tapeats  detritai  zircon 
ID-TIMS  data  from  ^Stewart,  et  al.  (2001 ),  and  two  LA-MC-ICP-MS  datasets  following 
the  procedure  of  Gehrels,  et  al.  (2008)  from  a  repeat  of  the  original  Stewart  sample 
and  a  second  Tapeats  Sandstone  sample  (Gehrels,  unpublished  data). 


Liverpool  Land,  east  Greenland  and  closely  associated  with  widespread 
migmatitic  textures  was  selected  for  a  first  application  of  the  small- 
volume  technique  to  a  suite  of  complicated  zircons.  With  a  primary 
phase  assemblage  of  l<yanite+ garnet +biotite+plagioclase+K-feldspar, 
the  sample  also  includes  abundant  rutile  and  accessory  monazite  and 
zircon.  Separated  zircons  are  primarily  30-80  pm  in  diameter  with  a 
rounded,  colorless  and  vitreous  crystal  habit  yielding  a  “beadlike” 
appearance;  <10%  have  yellow-tan  cores  and  are  occasionally  tabular 
and  up  to  100  pm  in  length.  CL  images  of  mounted  grains  display  a  range 
of  bright  and  dark  cores  exhibiting  oscillatory  through  isochemical 
growth  zoning,  whereas  rims  typically  <20  pm  in  width  display 
homogenous  CL  response  from  grain  to  grain. 

To  characterize  the  detritai  signature  and  the  timing  of  high- 
temperature  metamorphism  in  this  sample,  10  grains  were  selected 
randomly  for  U-Pb  age  mapping  by  the  small-volume  technique  with 
32  laser  pulses  (Electronic  Supplement,  Table  ES4A,  ES4B).  2a 
systematic  error  accumulated  over  the  course  of  the  sample  session 
was  1.4%,  and  10  analyses  on  secondary  zircon  standard  R33  measured 
during  the  analysis  session  yield  a  206p|^^238fj  424.7  ±7.2  Ma 


(total  2a  error  and  within  error  of  the  published  ID-TIMS  age)  and 
provide  a  quantitative  result  supporting  the  external  reproducibility  of 
the  analysis  session.  Working  systematically  from  the  rims  to  the  cores 
in  order  to  avoid  contaminating  possibly  younger  rims  with  ejection 
debris  derived  from  the  cores,  between  eight  and  14  analyses  were 
made  in  each  grain  making  sure  to  analyze  material  from  the  various 
chemical  domains  identified  in  CL.  Plotting  all  analyses  on  a  concordia 
diagram  yields  concordant  ages  spanning  the  Archean-Mesoproter- 
ozoic  with  discordant  analyses  defining  a  wedge-shaped  data  field 
converging  toward  a  lower  intercept  in  the  Silurian  (Pig.  12 A).  U/Th 
ratios  are  typically  <10  for  core  analyses,  whereas  rim  and  analyses 
with  ages  younger  than  -500  Ma  yield  higher  U/Th  ratios  up  to  126 
(Pig.  12A,  inset).  Calculating  upper-intercept  ages  after  Ludwig  (2003) 
from  analyses  within  individual  grains  yields  Model  1  ages  ranging 
from  1180-2840  Ma  with  age  errors  that  vary  from  1-11%  (95% 
confidence  +2a  systematic  error).  Under  the  assumption  that  all 
grains  grew  younger  zircon  rims  at  the  same  time,  analyses  from  all 
grains  with  206p|^y238y  207p|^y235y  <500  Ma  were  pooled  to 

calculate  a  lower-intercept  age.  206p|^y238y  207p|^^235y 

434.7  +10.4/-6.3  and  433.7 +  7.9/-6.8  Ma  (95%  confidence  plus  2a 
systematic  error),  respectively,  were  then  calculated  using  the  TuffZirc 
algorithm  of  Ludwig  (2003),  designed  to  identify  statistical  outliers 
and  calculate  ages  that  are  minimally  affected  by  either  inheritance  or 
Pb-loss  (Pig.  12B).  Nine  out  of  the  10  mapped  grains,  each  with  cores  of 
different  age,  yielded  at  least  one  analysis  that  was  accepted  by  Tuff¬ 
Zirc  and  included  in  the  age  calculation,  and  demonstrates  the  ability 
of  this  technique  to  target  specific  age  domains  in  intricately  zoned 
crystals.  Although  the  trimmed  dataset  suggests  a  minor  inherited 
component  on  concordia  plot  (Pig.  12B),  the  data  define  a  single 
population  and  further  trimming  of  the  dataset  is  not  warranted  given 
the  assigned  errors.  Archean-Mesoproterozoic  detritai  zircon  ages 
and  a  high-temperature  Silurian  metamorphic  event  are  consistent 
with  detritai  signatures  and  widespread  anatexis  from  440-425  Ma 
observed  in  the  Krummedal  Group  in  east  Greenland  farther  west  of 
Liverpool  Land  (Kalsbeek  et  al.,  2000;  Watt  et  al.,  2000;  Jones  and 
Strachan,  2000;  Kalsbeek  et  al.,  2001;  White  and  Hodges,  2003). 

5.  Conclusion 

The  small-volume  U-Pb  geochronology  technique  at  the  Arizona 
LaserChron  Center  microsamples  and  dates  zircon  volumes  as  small  as 
12-14  pm  in  diameter  by  4-5  pm  in  depth,  comparable  to  sample 
volumes  for  U-Pb  geochronology  by  secondary  ion  mass  spectrometry 
(Ireland  and  Williams,  2003).  Low  Pb  yields  produced  by  the  small  laser 
Spot  diameter  and  a  slow  laser  pulse  rate  require  using  Channeltron 
detectors  to  monitor  all  Pb  peaks  coupled  with  a  total  count  integration 
counting  method  to  calculate  U  and  Pb  peaks  and  isotope  ratios. 
Measurement  errors  are  calculated  as  a  function  of  the  number  of  counts 
and  observed  206p|^^238y  on  Sri  Lanka  standard  zircon.  Using  this 
technique,  precise  and  accurate  ages  on  eight  Mesoproterozoic- 
Cretaceous-aged  secondary  zircon  standards  were  reproduced  within 
2%  of  published  ID-TlMS  ages;  a  ninth  Miocene-aged  standard  produced 
an  age  - 1  Ma,  or  3-4%,  too  young.  In  two  initial  applications  of  this 
technique,  small-volume  U-Pb  geochronology  was  used  to  determine 
the  detritai  zircon  signature  of  the  Bright  Angel  Shale  and  to  make  U-Pb 
age  maps  of  individual  zircons  to  unravel  the  detrital-metamorphic 
history  of  a  granulite-facies  paragneiss  from  east  Greenland.  In  addition 
to  the  small  spot  diameter,  the  shallow  pit  depth  is  essential  to  the  utility 
of  the  small-volume  technique  because  it  ensures  that  the  sampled 
volume  corresponds  directly  to  chemical  zones  observed  in  CL  and 
minimizes  the  probability  that  multiple  chemical  zones  will  be  sampled 
at  greater  depths  within  the  sample.  The  ability  of  this  method  to 
overcome  difficulties  associated  with  acquiring  meaningful  U-Pb  ages 
on  complexly  zoned  zircons  by  targeting  single  age  domains  <15  pm, 
combined  with  the  rapid  rate  of  data  collection  (-60  analyses  per  hour) 
and  non-destructive  nature  of  the  technique,  should  prove  useful  for  a 


Please  cite  this  article  as:  Johnston,  S.,  et  al..  Small-volume  U-Pb  zircon  geochronology  by  laser  ablation-multicollector-ICP-MS,  Chemical 
Geology  (2008),  doi:10.1016/j.chemgeo.2008.11.004 


ARTICLE  IN  PRESS 


S.  Johnston  et  al  /  Chemical  Geology  xxx  (2008)  xxx-xxx 


11 


Fig.  12.  A)  Concordia  plot  and  U/Th  data  for  ail  analyses  acquired  over  the  course  of  making  U-Pb  age  maps  on  10  grains  from  the  Greenland  granulite-facies  pelite.  As  an  example  of 
how  upper  intercept  ages  are  calculated,  analyses  from  grain  6  are  shown  as  filled  ellipses  with  analysis  numbers  on  the  concordia  plot  that  correspond  to  spot  locations  on  the  inset 
CL  image.  B)  Concordia  plot  of  all  analyses  with  206p|^y238y  207p|^y235y  <5qq  jyi^  adjacent  to  Tnjflxrc  results;  discarded  analyses  are  indicated  with  open,  grey  ellipses  and  grey 

error  bars,  respectively. 


wide  variety  of  geologic  applications,  and  confirms  the  utility  of  LA-MC- 
ICP-MS  as  a  tool  for  high  spatial  resolution  U-Pb  zircon  geochronology. 
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